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DISLOCATION PATTERNING UNDER 
INTERNAL AND EXTERNAL STRESSES 

Mahesh Kumar Sahu 


SYNOPSIS 


Dislocations in one- and two- dimensions have been simulated by a computational 
technique similar to molecular dynamics. The model employs an array of straight 
and pure edge dislocations on parallel slip planes. Externally applied stresses and 
interaction forces dm 4 to other dislocations are integrated using a time stepping 
scheme. In the case of dislocations moving under internal stresses, it is found that 
the present model reproduces the analytical results accurately. A random two- 
dimensional array of dislocations found break up into two walls under the action of 
internal stresses. In addition, these walls become wavy under the action of a cyclic 
load. These results may be useful in understanding the formation of the persistence 
slip bands under fatigue loading. 
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Nomenclature 


u, ; Displacement in i direction, 
u : Resultant displacement. 

<7 tJ : Stress vector acting in i direction on j plane . 
t'u : Strain vector acting in i direction on j plane . 

G : Shear modulus. 

A : Lame constant. 
b : Burgers vector. 

b, : Burgers vector’s xth component of dislocation. 
v : Poisson’s ratio. 

(r , 0) -.Radial parameters of two-dimensional coordinate. 
F : Force acting on dislocation, 
r : Shear stress. 
v : Velocity of dislocation. 

A t : Time step of integration. 

a*, : ar-eoordinate of dislocation on step i. 

a l :Stress field due to x components of Burgers vector. 

er2 :Stress field due to y components of Burgers vector. 

T ex t • External shear stress. 

r m : Mean of cyclic stress. 

r 0 : Amplitude of cyclic stress. 

T : Periodicity of cyclic stress. 
x + ^-coordinate of positive dislocation. 
x” :ar-coordinate of negative dislocation. 
x mm : Half width of dislocation distribution. 



Chapter 1 
Introduction 


1.1 Introduction 

1.1.1 Existence of Dislocations 

Although there are many techniques now available for the direct observation of 
dislocations, the existence of these line defects was deduced by interference in the 
early stages of dislocation study. Strong evidence arose from attempts to reconcile 
theoretical arid experimental values of the applied shear stress required to plastically 
dv.form a single crystal. The striking difference between prediction and experimental 
was accounted for by the presence of dislocations independently by Orawan, Polanyi 
and Taylor in 1934. 

In recent years, it has been possible to produce crystals in the form of fibers of a 
small diameter, called whiskers, which have a very high degree of perfection. These 
witiskers are sometimes entirely free of dislocations and there strength is close to 
the theoretical strength. 

Dislocations are mainly divided in two types. 

• Edge dislocation 

• Screw dislocation 

Practically the dislocations are in mixed state of those two types. 
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1.1.2 Statement of The Problem 

The experimental studies of dislocation arrangements in deformed crystals are com- 
promised by the fact that most examinations must be performed ex post facto. 
Consequently such observations cannot address questions of dynamic dislocation 
interactions and arrangement as functions of dislocation content, character, thermo- 
mechanical conditions and crystal geometry. Since few materials are amenable to 
dynamic observations of dislocations during deformation, computational modelling 
must be employed to determine the nature and geometry of dislocation arrange- 
ments under these conditions. 

Modelling the evolution of steady-state dislocation arrangements under the action of 
internal and external stresses is a many-body problem. A complete three-dimensional 
formulation involves consideration of changes of dislocation shape both within and 
out of the slip plane as well as appropriate formulation of boundary conditions which 
do not introduce artifacts into the solutions. Such formulations are limited by the 
capacity of available computing resources to problems involving relatively few dis- 
locations. 

Simulation of dislocation arrangements in two dimensions substantially reduces the 
computational requirements of the problem while revealing many useful aspects 
of the evolution of dislocation structure and various thermomechanical conditions. 
This permits many more variables to be introduced into the problems without in- i 
earring an unacceptable penalty in computational requirements. The results of such i 
computations naturally must be interpreted with the simplifications in mind. i 

Dislocations in deformed metals tend to cluster into various kinds of dense regions j 
of high dislocation density separated by the dislocations of low dislocation density. !. 
Observation of dislocation structures in the electron microscope reveals that tend to 
tangle and to form cells. In worked metals of high stacking fault energies, dislocation 
are present in planar arrays of high energy, which constitutes a strong driving force j 
for recrystallization. J 

Low- and high-temperature deformation conditions are defined by the relative im- j 
portance of glide and climb. Dislocation arrangements resulting from low- and j 
high- temperature deformation exhibit significant differences. At high temperatures 
dislocations gain an additional degree of freedom through their ability to climb. The 
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rate of deformation of high-temperatures is primarily controlled by atomic diffusion, 
which is negligible at low temperatures. After removal of applied stress, dislocations 
relax into configurations determined by the initial configurations, the relative mobil- 
ity of glide and climb(We have neglected in our code) and the mutual forces among 
dislocations. 

The present study simulates the dislocations both in one-dimensional and two- 
dimensional. low-temperature condition , where only glide motion is permitted. Effects 
of cyclic stresses are also studied. 

1.2 Dislocation Theory and Collective Effects 

1.2.1 Main Object and Achievement 

The main object of dislocation theory is the prediction of the mechanical properties 
of bulk materials, or at least their understanding, in terms of the physical properties 
of the material considered and of its structural defects. To date this remains a dis- 
tant aim despite many detailed studies performed in carefully controlled conditions, 
in particular on single crystals. 

A major achievement of dislocation theory is concerned with the interaction of a mo- 
bile dislocation segment with a localized defect which is assumed to be distributed 
either periodically or at random.The plastic strain rate is obtained by summing up 
over the total no. of mobile segments the individual properties for cutting through 
such obstacles, which is equivalent to assume that all events are independent. This 
leads within a good numerical accuracy to the prediction of the thermal component 
of the yield stress. 


1,2.2 Difference in Local and Average Properties j 

I 

Models which assume plastic flow to be homogeneous assimilate the local strain rate; 
to the average strain rate and, therefore, identify two quantities which differ by one! 
to two orders of magnitude, at least at moderate strains. ! 

More generally, strain nonuniformities may occur at various scales of observations,! 
macroscopic, mesoscopic(i.e., optical or metallographic) and microscopic.Structural 
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instabilities can in principle be investigated in terms of dislocation theory, while 
geometrical instabilities (for instance necking or several types of shear banding) are 
better described within the framework of micromechanics. 

In both eases very basic problems such as the conditions for the occurrence of strain 
nonuniformities, their scaling laws and the kinetics of their further development are 
still not well understood. 


1.2.3 Dislocation patterning 

At moderate and large strains an increasing density of immobile or slow disloca- 
tions get stored in the deforming crystals, leading to an increased glide resistance, 
possibly to long range internal stress and strain hardening. As revealed by transmis- 
sion electron microscopy these dislocations from various of bi- or three-dimensional 
arrangements, e.g. wall and channel structures, where dislocation-rich regions al- 
ternate more or less regularly with dislocation-poor regions (Louchet and Brechet, 
1988; Amodeo, 1988; Hansen and Kuhlmann-Wilsdorf, 1986) 

The relation between microstructure and the strain hardening properties of solids 
has been the object of many phenomenological models. In the simplest ones the dis- 
location densities are considered and their evolution with time or strain is through 
a balance between production and storage or annihilation events. More elaborated 
models assume a particular type of spatial distribution for the mobile and stored 
dislocations, and consider it as the smallest elementary volume whose mechanical 
properties are identical to those of the bulk crystal. 

The underlying questions, which are still at debate are as follows: 

* What are the conditions of occurrence for such type of ordered or semi-ordered 
microstructures, which we further refer as “pattern”? 

* What are the material parameters governing their geometry, their evolution, 
their transformation during straining ? 

* How they may be related to strain hardening properties? 
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• Is the formation of a field of long range internal stresses necessarily associated 
with patterning? 

1.3 Progress in Dislocation Study 

Although this brief description of dislocation theory and of its limitation may seem a 
bit schematic, the need for a consistent modeling of collective dislocation effects and 
of their consequences regarding strain hardening and strain localizations has been 
recognized very early. The tools, both conceptual and numerical, needed to develop 
such an approach have only been made available in the last two decades through 
the investigation of the dynamic behavior of ensembles of interacting populations. 

The application to dislocation theory was initiated in the mid-eighties and it has 
yielded a few promising primary results. Only a few idealized situations have been 
investigated to date and this domain is still in its infancy. The connection between 
collective dislocation effects and patterning phenomena occurring in other systems 
brought far from thermodynamic equilibrium is, however , firmly established. It 
can he illustrated by recent aehievement,s(Brechet, 1987; Lepinoux, 1987; Amodeo, 
1988; Schiler, 1989;) j 

! 

1.4 Research Objectives j 

The field of dislocation is relatively recent topic for scientists.This stream is opened 
for many innovations. However numerous studies have been conducted of dislocation | 
arrangements formed under thermomechanical conditions. The dynamic behavior of j 
individual dislocations has been studied theoretically, by computer simulation and j 
experimentally. Some studies based upon automata and molecular dynamics simu- j 
lations for many dislocations are available. However, these model do not consider j 
the true long-range nature of the dislocation-dislocation interactions, which control | 
many of the features of dislocation microstructures. Also these models restricted to | 

single slip system. j 

The objects of the research include to understand the characteristics of dislocations 
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when they are simulated on the basis of conventional ways, like stress field of other 
dislocations. 

1. What will be the final arrangements of dislocations when they will be simulated 
in resultant stress field of other dislocations? 

2. The effects of external stresses on dislocations configurations . 

The above 1 and 2 steps can be studied for one- as well as two-dimensional array of 
dislocations. 

1.5 Overview of the Study 

This thesis consists of total six chapters. First chapter includes the primarily back- 
ground of the dislocation, short introduction and objects of this research. Chapter 
two is containing the very important part: formulation of stress field of edge and 
screw dislocations. Computational procedure and flow chart is also elaborated in this 
chapter. In chapter three the one-dimensional simulation is explained with compu- 
tational and analytical method.The internal stress and externally applied stresses 
are include*! for computational steps.Chapter four is similar to chapter three except 
that this explains the two-dimensional simulation.In chapter five the all important 
results are included and the graphs are explained.In chapter six results are concluded 
anti scope for future work is pointed out. 

1.6 Figures of Dislocations 

Many types of dislocations are shown on next two pages: 







Mixed dislocation: 

(a)The curved dislocation SME is pure edge at E and pure screw at S. 

(W Burgers vector b of dislocation XY is resolved into pure edge componen t bl 
and pure screw component b2 


Figure 1.2: Mixed Dislocation 




Chapter 2 

Dislocation Dynamics 


2.1 Stress Field of Dislocations 

2.1.1 Introduction 

The atoms in a crystal containing a dislocation are displaced from their perfect 
lattice sites, and the resulting distortion produces a stress field in the crystal around 
the dislocation. The dislocation is therefore a source of internal stress in the 
crystal. The stresses and strains in the bulk of the crystal are sufficiently small 
for conventional elasticity theory to be applied them. This approach only to be > 

ceases valid at positions very close to the center of the dislocation. Although most 

i 

crystalline solids are elastically anisotropic, it is much simpler to use isotropic ; 
elasticity theory. This still results in a good approximation in most cases. j 

From a knowledge of the elastic field, the following important characteristics can be j 
obtained: j 

j 

• The energy of the dislocation. ! 

• The force it exerts on other dislocations. 

! 

• It’s energy of interaction with point defects. j 

The elastic field produced by a dislocation is not affected by the application of stress 
from external sources: the total stress on an element of the body is the superposition ! 

of the internal and external stresses. 
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2.1.2 Elements of Elasticity Theory 

The displacement of a point in a strained body from its position in the unstrained 
state is represented by the vector, 

U = [«*, Uy, u z ] (2.1) 



Displacement of F to P’ by displacement vector u 
Figure 2.1: Displacement of a point 


The components u x , u„, u t represents projections of u on the x, y,z axes, as 
shown in Fig. 2.1. In linear elasticity , the nine components of strains are given in 
terms of the displacement components thus: 

•Normal strains are: 


dtlr $Uy 

m w eyv = w 


e xz = 


du z 

dz 


(2.2) 


•Shear strains are: 
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(a) 



(b) 


(a)Pure shear and (b)simple shear of an area element in the xy plane. 

Figure 2.2: Pure Shear 



Normal strains represent the fractional change in length of elements parallel to 
the x, y and z axes respectively. The shear strains can be understood by the Fig. 
2.2(a), in which the small area element ABCD in the xy plane has been strained 
to the shape AB'C'D' without change of area.The angle between side AB and AD 
initially parallel to x and y respectively has decreased by 2e*j,. By rotating, but not 
deforming, the element as in Fig 2.2(b), it is seen that the element has undergone 
a simple shear. The simple shear often used in engineering practice is 2e*j„ as 
icatcd. 

actional change in volume A, known as the dilation , is therefore: 

A s= = (e xx 4- e yy + e„) (2.6) 

* independent of the orientation of the axes x,y,z. 
relationship between stress and strain in linear elasticity is Hooke’s Law in 
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which each stress component is linearly proportional to each strain. For isotropic 
solids, only two proportionality constants are required: 

&XX ^GGxX ~t* "1" Gyy Ggz ) 

®yy 2£?eyy + A (e xx + 6yy -)- e zz ) 

CT zz H” X{g X x d™ C yy d™ G zz ^ 

°*U 2Gv-xy n ®ys = <7 ~ x = 2Ge zx (2-7) 

A and G are the Lame Constants, but G is more commonly known as the shear 
modulus. 

2.1.3 Stress Field of Straight Dislocation 

Screw Dislocation: 

The elast ic Held around an infinitely-Iong, straight dislocation can be represented 
in terms of a cylinder of elastic material. A radial slit LMNO was cut in the cylin- 
der parallel to the 2 -axis and the free surfaces displaced rigidly with respect to each 
other by the distance 6, the magnitude of the Burgers vector of the screw dislocation, 
in the 2 -direction. 

'Hie elastic Held in the dislocated cylinder can be found by direct inspection. 
First, it is noted that there are no displacements in the x and y directions: 

u x = u y = o (2.8) 

Secondly, the displacement in the z direction increases uniformly from zero to b as 
0 increases from 0 to 2tt: 

It is then found from equations (2.2) to (2.5) that 

e xx = €yy s® &zz = &xy = &yx ^ O 


(2.9) 
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z 

A 




(a) Screw dislocation AB formed in a crystal. 

(b) Elastic distortion of a cylidrical ring simulating 

the distortion produced by the screw dislocation AB. 

Figure 2.3: Screw Dislocation 


_ b y _ b sin 6 

( IZ < 4tt (x 2 + y 2 ) Arc r 

- — ^ - — b cos ^ 

e ' vz ~~ 4jt (x 2 + y 2 ) ~ An r 
From equations (2.7) and (2.10), the components of stresses are 

Oxx — O'yy — <J iZ = <?xy — C 7y X — 0 

_ Gb y _ Gbsin$ 

Oxt - = - 2^ ^2 + y2 ) - - 25? r 

~ ^ cos ^ 

<7y * ” ” 2?r (a 2 + y 2 ) 2n r 


( 2 . 10 ) 


( 2 . 11 ) j 
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We can conclude from stress and strain equations ( 2.10) and (2.11) that: 

• The elastic distortion contains no tensile or compressive components and con- 
sists of pure shear. 

• The stresses and strains are proportional to 1/r and therefore diverge to infin- 
ity as r -4 0. Real crystals are not hollow as taken in this analysis, and so as 
the center of a dislocation is approached, elasticity theory ceases to be valid 
and a nnn-linmr atomistic model must be used. The region within which the 
linear elastic solution breaks down is called the core of the dislocation. 

Edge Dislocation: 




ho 

(a)Edge dislocation formed in a crystal . 
tb) Elastic distortion of a cylindrical ring simulting the 
distortion produced by the edge dislocation in (a) 

Figure 2.4: Edge Dislocation 

The stress field is more complex than that of a screw but can be represented 
in an isotropic cylinder in a similar way. Considering the edge dislocation in Fig. 
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2.4(a) elastic strain field can be produced in the cylinder by a displacement of the 
faces of the slit by a distance b in the a>direction(Fig. 2.4(b)) The stresses are: 


Vxx = —Dy 


(3a; 2 + y 2 ) 
(a; 2 + y 2 ) 2 


a yy 


= Dy 


(a ; 2 - y 2 ) 

(x 2 + y 2 ) 2 


&xy 


&yx 


= Dx 


(a; 2 - y 2 ) 
(a; 2 + y 2 ) 2 


Ozz = v{o xx + <Jyy) 

Cjcj SB 0 ZX = CTy z == CT z y = 0 


( 2 . 12 ) 


where 

D = — — — 

27 r(l - u) 

We can say from above equations that: 

• The stress field lots, therefore, both dilational as well as shear components. 
The largest, normal stress is a xx which acts parallel to the slip vector. 

* The elastic solution has an inverse dependence on distance from the line axis 
ami breaks down when x and y tend to zero. It is valid only outside a core of 
radius r&. 

2.2 Forces on Dislocations 

When a sufficiently high stress is applied to a crystal containing dislocations, the 
dislocation move and produce plastic deformation either by slip or at sufficiently 
high temperatures by climb. The load producing the applied stress therefore does 
work on the crystal when a dislocation moves, and so the dislocation responds to 
the stress as though it experiences a force equal to the work done divided by the 
distance it moves. Only glide force is considered here, climb force is neglected. 
Consider a dislocation moving in a slip plane under the influence of a uniform 
resolved shear stress r (Fig. 2.5). When an element d l of the dislocation line of 
Burgers vector b moves forward a distance ds the crystal plane above and below 
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The displacement ds of an element dl in its glide plane 

Figure 2.5: Force on Dislocation 


the slip plane will be displaced relative to each other by b. The average shear 
displacement of the crystal surface produced by the glide of dl is: 

(2.13) 

Where .*1 is t he area of the slip plane. The external force on this area is At , so the 
work done when the element of the slip occurs is 


dW = Ar(^~~)b (2.14) 

the glide farce on a unit length of dislocation is defined as the work done when unit 
length moves unit distance. Therefore, 


F = 


dW dW 
dsdl = d^4 


(2.15) 


The stress t is the shear stress in the glide plane resolved in the direction of b and 
the glide force F acts normal to the dislocation at every point along its length. 
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2.3 Velocity of Dislocation 

Dislocations move, by glide at velocities which depend on the applied shear stress, 
purity of crystal, temperature and type of dislocation. A crystal containing freshly 
introduced dislocations usually produced by lightly deforming the surface, is sub- 
jected to constant stress pulse for a given time. From the positions of the dislocations 
before and after stress pulse, the distance each dislocation has moved, and hence 
the average dislocation velocity can be determined. 

Experimentally the dislocation velocity was measured and observed that dislocation 
velocity was a very sensitive function of the resolved shear stress.In the range of 
velocities between 10“ 7 and 10 -1 cms -1 , the logarithm of the velocity varies linearly 
with the logarithm of the applied shear stress, mathematically: 



where v is the velocity, r is the applied shear stress resolved in the slip plane, r 0 is 
the shear stress for v = 1 cms~ l and n is a constant and was found experimentally 
to be ~ 25 for lithium fluoride.lt must be emphasized that equation (2.16), is purely 
empirical and implies no physical interpretation of the mechanism of dislocation 
motion. 

The stress dependence of dislocation velocities varies significantly from one mate- 
rial to another. For a given material, the velocity of transverse shear wave propaga- 
tion is the limiting velocity for uniform dislocation motion. However damping forces 
increasingly oppose motion when the velocity increases above about 10" 3 , and thus 
n in equation (2.16) decreases rapidly in this range. 

Studies on face centered cubic and hexagonal close packed crystals has shown that 
at the crystal resolved shear stress for macroscopic slip, dislocation velocity is ap- 
proximately 1ms" 1 (10“ 3 of a shear wave velocity) and satisfies 

v oc Ar m (2.17) 

Where A is the material constant and m is approximately 1 at 300 K in pure crystals. 
We want the relation between velocity of dislocation and force acting on dislocation 
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for simulation purpose. So we will replace r in terms of F, 


v (x AF/b 


v = 



(2.18) 


Where C is a proportionality constant. Now we can put all constant terms to unit 
so that we can get convenient form for simulation, we can write the above equation 
as: 


v = F 


dx 
d t 


= F 


(2.19) 


2.4 Time Integration Scheme 

For getting the positions of dislocations on different time steps, difference scheme 
of numerical method can be used. For getting the position after first time step we 
have to use forward difference method. 
equation (2.16) can be written in the form of difference: 

2>t+l __ rp' 

At “ *’ 

Xi+i =Xi + A tFi (2.20) 

Where x i+ i and are the positions of dislocations at (i + l)th and ith time steps 
respectively, At is time-step and Ft is the force acting on dislocation at ith time- 
step. 

Now we can use for positions of further time steps with the central difference method 
as that is more accurate than forward difference method. 

Xj+l ~ -Ej-l p, 


t, , , = T.- , 4- 2A tFi 


( 2 . 21 ) 
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Graphical representation of forward difference method. 

Figure 2.6: Difference Method 


The equation (2.17) and (2.18) is applicable when At is between 0 and 1 and will 
converge when At tends to 0. 

2.5 Forces between Dislocations 

The basis of the method used to obtain the force between two dislocations is the 
determination of the additional work done in introducing the second dislocation into 
the crystal which already contains the first. Consider two dislocations lying parallel 
to the z-axis. (Fig 2.7). 

The total energy of the system consists of: 

(a) the self energy of the dislocation I. 

(b) the self energy of dislocation II, and 

(c) the elastic interaction energy between I and II. 

The interaction energy Eint is the work done in displacing the faces of the cut which 
creates U in the presence of stress field I. The displacements across the cut are 6*, 
by, b z the components of the Burgers vector b of II. By visualizing the cut parallel 
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Figure 2.7: Interaction between two edge Dislocations 
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to the x or y axis, two alternative expressions For E int per unit length of II are 
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ft 
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+ 6^cr^)dar 


r oo 

E { „, = / (4 X^XX + b y a yx + 6 z (7 zx )dy 

Jy 


( 2 . 22 ) 

(2.23) 


where the stress components are those due to I. The interaction force on II is ob- 
tained simply by differentiation of these expressions, i.e. F x - -dE int /dx and 
F y = - dEi nt /dy . 

For the two parallel edge dislocations with parallel burger vectors shown in Fig. 2.7, 
by ss b,c ss 0 and b x s = 6, and the components of the force per unit length acting on 
II are therefore 

F x = o xy b, Fy « a xx b (2.24) 

where a xy and a xx are the stresses of I evaluated at position (a:, y) of U. The forces 
if tt is a negative edge. Equal and opposite forces act on I. F x is the force 


i 


! 
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in the glide direction and F y the force perpendicular to the glide plane. Substituting 
from equation (2.12) gives 


Gfr 2 x{ s 2 - y 2 ) 
27r(l — v) (x 2 + y 2 ) 2 


(2.25) 


G6 2 y( 3x 2 + y 2 ) 
y 27r(l — z/) (a: 2 + y 2 ) 2 


(2.26) 


Since edge dislocations can move by slip only in the plane contained by the disloca- 
tion line and its Burgers vector, the component of the force which is most important 
in determining the behavior of the dislocations in Fig. 2. 7 is F x . 

We can conclude from equations (2.25) and (2.26) that: 



Figure 2.8: Force between parallel edge Dislocations with parallel Burgers vectors. 


• for x > 0,F X is negative (attractive) when x < y for dislocations of same sign 

* - .ff. — /vf sifrn. 



2.6 Unit Dislocation Field 


22 


• for x < 0 ,F X is positive (attractive) when x > — y for same sign and x < —y 
for the opposite sign. 

F x is plotted against x, expressed in units of y , in Fig.2.8. It is 0 whenx = 0, ±y, ±oo, 
but of these, the positions of stable equilibrium are seen to be x = 0, ±oo for edges 
of the same sign and ±y if they have the opposite sign. 

It follows that an array of edge dislocations of the same sign is most stable when 
the dislocations lie vertically above one another as in Fig.2.9(a). Furthermore, edge 
dislocations of the opposite sign gliding past each other on parallel slip planes tend 
to form stable dipole pairs as in Fig. 2.9(b) at low applied stresses. 



Figure 2.9: Stable positions for two edge dislocations of (a)the same sign and (b) 
opposite sign. 


Comparison of the glide force F x in equation(2.15) shows that since a xy is the shear 
stress in the glide plane of the dislocation II acting in the direction of its Burgers 
vector, equation (2.15) holds for both external and internal sources of stress. 

2.6 Unit Dislocation Field 

Let a dislocation is placed at (a:, y) coordinate and we want the stress field at origin. 
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in x direction (6 r ), is. 


-i _ ~y ( 3a;2 + y 2 ) 
^ ( X 2 + y 2)2 

_ y(x 2 - y 2 ) 

yy (: x 2 + y 2 ) 2 
, _ x(z 2 - y 2 ) 

XJ/ ~ (x 2 + y 2 ) 2 


In terms of radial coordinates the above equation can be written as: 

sin 0(cos 26 + 2) 


ah 


®^-yy ~ 


a 1 xy — 


sin 6 cos 26 


cos 6 cos 26 


From Fig. 2.10, we can got, the stress field for b y by replacing in x 
So the stress field will be. 

-x( 3y 2 + x 2 ) 


<t2, 


yy ~ (x 2 + y 2 ) 2 
x(x 2 - y 2 ) 


a2 xx 
a 2„ 


( x 2 + y 2 ) 2 
y(x 2 - y 2 ) 


■' ! ' x ~ (x 2 + y 2 ) 2 
In radial coordinates and in form of a 1 


a 2 


yy 


cos #(2,0 - cos2fl) 

r 

a2 xx = a 1 xy 

02yx a lyy 

So the resultant stress field by b (b x , b y ): 

a xx ~ b x a I®* + bya2 xx 


(2.27) 


(2.28) 
y and y -4 x. 


(2.29) 


(2.30) 
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(a) 

By rotating (a): 


c r 
x 


X 


T 


(b) 

Which is similar to Fig(b) . 

If we replace if X=-Y and Y=X 


Figure 2.10: Stress Field for Different Orientation. 


(Tyy — b X (Tlyy + by(X2yy j 

&xy = b x <yl X y + by&2 X y (2.31) 

Let us assume that there are n dislocations for simulation. Each dislocation is j 
influenced by the resultant stress field of all other dislocations. So the expression j 
for force acting on ith dislocation will be (from equation 2,24): 

n 

Fi = bi %(*» i) when j # i (2.32) 

i = i 

where r xy (i,j) is the stress field of dislocation j on position of ith dislocation, 6, is 
Burgers vector of ith dislocation and Fi is the force acting on ith dislocation in x 
direction due to stress field of all other j dislocations. 
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2.7 Computational Procedure 

The specific computational procedures performed by dislocation simulation code are 
explained in following steps: 


1. input dislocation parameters like coordinates, Burgers vectors ,time and 
time-step for simulation. In case of applied stress, parameters of applied 
stress like periodicity, mean value and amplitude value of stress are also 
given as input. 

2. The total force applied on a single i dislocation by all other j dislocations is 
calculated. 

• Firstly the distance between dislocations in form of x and y components 
is calculated and that is got finally in the form of radial parameters 

M). 

• From equations (2.28) to (2.31) (from radial distance parameters (r, 9)) 
tin* stress field on position of dislocation i is calculated.In case of 
externally applied shear stress, this stress is simply added to internal 
stress field. 

• By multiplying this stress field with Burgers vector we calculate The 
force (From Equation 2.32) experienced by the ith dislocation.By doing 
this process in a loop and adding the calculated force in previous value 
of force ,wc calculate the total force applied by all other dislocations. 

• From equations 2.20 and 2.21, we can get the new coordinate position of 
dislocation firstly by forward difference, for further steps we will use 
central difference method for getting more accuracy in calculation. 

• These calculation are done for every dislocation and we get a new set of 
coordinates of dislocation. 

3. These steps are iterated with very small time-step and we get the new 
positions of dislocations successively with further time steps. 
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These data in form of coordinates of dislocations is checked with data obtained 
analytically. We observed that the computational data is converging properly with 
analytical data when time-step is taken very less (in the order of~ 0.001). 

2.8 Flow Chart for Time Integration Method 



/ INPUT AS DISLOCATION 
PARAMETERS, TIME PARA- 
METERS.APPLIED STRESS 
PARAMETERS IF NEEDED 



CONTmm ON NEXT PAGE 
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Figure 2.11: Computational Flow Chart 







Chapter 3 

One Dimensional Simulation 


It should ho noted that in this thesis we are treating only unit pure edge 
dislocations , however the dislocations are of both positive and negative type. 


3.1 Case I:-Only with Internal Stress 

For evaluating the behavior of two dislocations, we will put positive unit 
dislocation at .r a distance while a negative unit dislocation is placed on the -xq 
distance from origin on x-axis, as shown in Fig.3.1. 


3.1.1 Analytical Solution 

Here only two dislocations are in a line for interacting with each other. The 
positive dislocation will be influenced by stress field of the negative dislocation and 
vice versa. The st ress field of negative dislocation (from equation 2.27) will be 

1 


1 (x + - x~) 

with stress field keeping all constant terms unit. 

X+ and x~ shows x coordinate of positive dislocation and negative dislocation 
respectively. On putting s+ = x and x~ - -x , above equation can be written 

dx 1_ 

dt 2x 

On integrating: 


(3.1) 


(3.2) 


/ 


dx2x 


■/ 


dt 
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Figure 3.1: Two Dislocations in a line. 

,r 2 = -t + C Where C is integration constant (3.3) 

Boundary condition for positive dislocation is: 

For fss 0 x- x„ 

On applying this boundary condition to the equation 3.2 , we get the analytical 

solution in the following form. 

x = \Jxl - t (3.4) 

We have taken t.q = 1 , so the analytical solution is: 

s = (3.5) 


3.1.2 Computational Solution 

For iterating, the equation 3.2 can be written from difference scheme as following: 

ffj+l ~ 3»- 1 _ l 

2 X{ 

1 

%i+i = 1 + At— 


(3.6) 
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The Error in analytical and computational result is 

Error = x i+ \ - \jxQ-t (3- 7 ) 

;r, + i is the position of positive dislocation at time t 

Further, we will see that the computational results are converging with analytical 
results. 


3.2 Case H:-Internal Stress and Externally 
Applied Stress 


We will simply add the applied stress to the internal stress of dislocations. We are 
applying a sinusoidal stress, which varies as shown in Fig.3.2. 



Figure 3.2: External Shear Stress. 


The external stress can be written as, 


. 27rt 

Text = fm d" 7a p. 


Where t m , r„ and T are mean, amplitude and periodicity of the external shear 


stress. 
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3.2.1 Analytical Result 

Considering only external stress 


On integrating. 


dx(t) 

IT = ^ 

x(t) =C + T m t + ^(l- COS 


Where C is t ht* integration constant. 
Boundary condition: t = 0, x(t) = xq 


... . Tr a f 2i xt\ 

X{t) = Xq + T m t + — f i I “ COS — ) 

Including the internal dress, The total force will be: 

1 , , /_ , _ 2 nt\ 

dt ~ (x'-xi) bi \ m a T ) 

For positive dislocation's coordinate we can write above equation as, 


d:r + 

dt 

For negative dislocation, 

da:” 

dt 


1 


(x+ - a: - ) 


+ ( r TO + r a sin 


27TI 


1 


( 


T m + r a sin 


2rt' 


(x~ - or 1 -) 

By adding and subtracting equations 3.10 and 3.11 and putting: 

x* + 


and 


x + - x~ 


We get two equations in terms of £ and rj, 

2-S = — — + 2 T ex t 
dt £ 

2^2 = 0 

dt 

Boundary conditions: when t~0, ( — xo and 17 = 0 


(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 
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The equation 3.12 can be written as, 

J = ___ + Tea;t (3. 14 ) 

Above equation is a non-linear differential equation ,so we will use following 
iteration procedure for numerical solution. 


2&t 


1 

-xr- +r a sm 


2ttA t 
T 


When r m is taken 0. 

?„ . = f,- , + 2A« ( - A + t. sin (3.15) 

From above e<juation we can get successive values of £, which is equal ^-coordinate 
of positive dislocation. 


3.2.2 Computational Result 


The equat ion for computational result: 


:r tf i - .t,-i 

2 At 


= k 


1 . 27 Tt\ 

- + r m + r a s.n -jrj 


keeping r m U 

JVf i = *i~ i + 2A tbi ^ ~ + r a sin (3.16) 

Now we can get the plots with different set of values of r a and T for two 
dislocations. The graphs of externally applied stress condition for two dislocations 
are depicted in chapter-5. 


3.2.3 For many dislocations 

The dislocations of both positive and negative signs are arranged on x-axis as 
shown in Fig.3.3. Dislocations will interact with all other dislocations by the force 
in x-direetion only. In Internal stress case, The effective force between them will 
be only due to resultant r xy of all dislocations. 

The effective formula for many dislocations simulation is 

3 

x i+ i = Xi-i + 2Atk ^ r xy (i,j) 

i=l 


(3.17) 
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When j ^ i 



Chapter 4 

Two Dimensional Simulation 


4.1 Internal Stress 

4.1.1 Dislocation Patterning 

Firstly wo will discuss about the behavior of dislocations due to internal stress field 
only. The dislocations are arranged in two-dimensional field randomly as shown in 
Fig. 4.1. 

We are going to neglect the inter-slip motion, so we will neglect the force acting in 
it direction. The effective stress field will be only r xy . The dislocations will glide 
only in their slip planes. 

The force acting on any ith dislocation can be calculated from equations (2.31) 
and (2.32). 

bx(j ')° r l*y ^y(j) cr ^xy < j (4.1) 

Where j ^ i 

Where b x and b v are Burgers vector’s components of jth dislocation-al^ and cr2*j, 
are stress field because of b x and b v of jth dislocation in position of ith dislocation. 
The expression for al xv and cr2 xy can be got from equations (2.28) and (2.30). 

By differential scheme we can get the successive coordinates of dislocations in 
further time steps. By iterating the above procedure for time of integration, we get 
the data as the coordinates of all dislocations with different time steps. 
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Figure 4.1: Dislocations in two dimensional-field. 


The plot of dislocation-patterning with internal stress are given in chapter-5 

4.1.2 Dynamics of Dislocation Wall 

Later we will see that due to internal stress the dislocations form the walls of 
dislocations (vertical dislocation-patterning). Dislocation-wall move with time 
away from each other. The velocity of these both walls can be calculated by the 
following steps. 

1. Firstly we divided the whole area (on which dislocations are moving) in 
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many small n vertical strips. That can be written in mathematical form as: 

/\ T max min 

n 

2. Now we can count the no. of positive and negative dislocations in each Ax 
region. The region on which no. of dislocations is highest contains the 
dislocation wall. 

H. In different time step we can get the position of dislocation by step 2. So we 
can get the plot between x-position and time for walls of positive dislocations 
and negative dislocations. So we can estimate the velocity of dislocation 
walls. 


4.2 Internal and External Stress 


For including t he external strews we have to add the external stress term in the 
equation 2.31, equation will turn out to be. 


<lx 

dJ 



( 4 . 2 ) 


In terms of parameters of cyclic external stress like mean, amplitude and 
periodicity of the stress the above equation can be written as. 


dx 


= 6 , 


E / * *\ 

r X y(t,j) + 


»= 1 



+ r a sin 


2irt" 


Rest procedure is same as internal stress computational procedure. 


( 4 . 3 ) 


4.3 Computational Flow Chart 

Computational flow chart for velocity of dislocation-wall is drawn on the next page 

as Fig.4.2. 



4.3 Computational Flow Chart 
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Figure 4.2: Flow Chart For Velocity of Dislocation-Wall 







Chapter 5 

JL 

Results and Discussion 


The dislocations arc studied for their movement in their slip planes in influence of 
two types of stress fields as described earlier. 

1 . ( )n!v by internal stresses. 

2. By interna! as well as external stresses. 

We will see further tiie results for one-dimensional and two-dimensional array of 
dislocations. 

5.1 Simulation for Two Dislocations 

5.1.1 Only Internal Stress 

The dislocations are placed on x-axis.Dislocation-I will interact with other by 
stress field of disloeation-II and vice-versa. The result is shown in Fig. 5.1. The 
figure is showing computational as well as analytical results, which are overlapping. 
EXPLAIN ATION: 

1. The dislocations arc gliding towards each-other with time before annihilating. 

2. The time taken for annihilation depends on their initial positions. 
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5-1.2 Internal Stress and External Stress 

Wo want to know the behavior of dislocation-gliding in presence of external 
stresses. We will see the effect in gliding with different amplitude and periodicity 
of external stress.The mean stress is kept 0 in all results. These results are 
depicted from Fig. 5.2a to 5.2f 
EXPLANATION: 

1. The dislocations gliding is unaffected by very small value of cyclic stress. 

(Fig. 5.2a) 

2. When wo increase the amplitude of external stresses, It starts to affect the 
dislocation motion. The stress tends to drive the dislocation gliding. (Fig 

5.2b) 

3. When the external stress is very high.The dislocation prominently drives the 
dislocation motion.(Fig. 5.2d and 5.2f) 

5.2 Many Dislocations in a line 

When many dislocations that randomly distributed in a line are subjected to 
internal stress meet of the dislocations remain around their initial coordinates. 
EXPLANATION: 

1. the dislocations glide only in z-axis. 

2. One negative dislocation moves to +oo while one negative dislocation to -oo 
.( From Fig. 5.3a to 5.3c) 

3. Most of dislocations vibrate around their original coordinates. 

5.3 Two Dimensional Dislocation 
5.3.1 With Internal Stress Only 
EXPLANATION: 
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1. the dislocations are moving in their glide planes only. 

2. The dislocations are making mainly a two vertical clusters (walls) of positive 
dislocations and negative dislocations. Fig. 5.4c 

3. The dislocation walls are moving to opposite direction with time. 

4. From plot 5.6a we see that the velocity is decreasing firstly and later it is 
achieving approximate constant value. 

5. velocity of positive and negative dislocation-walls are showing similarity. (Fig. 
5.6 b) 

6. When ,r min is increased successively the straight vertical wall (pattern) get 
defused. This is clear from Fig. 5.4c, 5.7a and 5.7b, in which plots only xmin 
is increasing 

7. If there is intrinsic length scale then we see that the walls splitting into the 
form of a sin wave.If we increase the size of the intrinsic length scale 
successively, the wavelengh of the wave decreases (from Fig. 5.8a and 5.8b). 

5.3.2 With Both Internal- and External-Stress Field 

The related plots are from Fig. 5.9a to 5.9b. 

EXPLAIN ATiON: 

1, In ease of both stress fields the dislocations are getting their vertical pattern 
but with time they are being drived by external stress field . 

2. Dislocation-walls are unstable with time . 

5.4 Important Graphs and Plots 

The graphs and plots are on next pages: 


Figure 8.1: Positive dislocation’s gliding with time by internal stress only 
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Figure 5.2a: Positive dislocation’s gliding with external stress 





I 


I 


.4 Important, Graphs and Plots 


44 


*0 



I Cn 



Figure 5.2e: Positive dislocation’s gliding with external stress 




Figure 5.2f: Positive dislocation’s gliding with external stress 





Figure 5.3b: Many dislocations on x-axis under internal stress only 





Figure 5.3e: Many dislocations on x-axis under internal stress only 
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Figure 5.4a: 250 dislocations clustering with time 
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Figure 5.4b: 250 dislocations clustering with time 




Figure 5.4e: 250 dislocations clustering with time 
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Figure 5.5: Density of 250 Dislocations After Integration 




Figure 5.6a: Velocity of positive Dislocation Wall with 250 Dislocations 





Figure 5.6b: Velocity of Both Walls 
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Figure 5.7a: Many Dislocations clustering with time: with different xmin 
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Figure 5.7b: Many Dislocations clustering with time: with different xmin 
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Figure 5.8a: Many Dislocations clustering with time: with diffrent intrinsic lengh 
scale 
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Figure 5.8b: Many Dislocations clustering with time: with different intrinsic lengh 
scale 
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Figure 5.9b; Dislocation Distribution:External Stress 



Figure 5.9c: Dislocation Distribution:External Stress 



Figure 5.9d: Dislocation Distribution:External Stress 




Chapter 6 
Conclusion 


6.1 Conclusion 

In the present study, the dislocations for one dimensional and two dimensional are 
being studied in case of applying external shear stresses. The unit edge 
dislocations are being studied for only their gliding in their slip planes.the 
distribution of dislocations by internal stress finally, is very surprising in case of 
two-dimensional distributed dislocations. 

Some important conclusions may be pointed out as following: 

1. In case of two dislocations in a line the dislocations are attracting towards 
each other till they annihilate . 

• The dislocation glide is prominently influenced by the external stresses 
,the gliding of dislocation is vibrated on £-axis by periodicity of external 
stress, if the amplitude and periodicity of cyclic stress is high. 

• When the external stress is of low periodicity and low amplitude the 
vibration is very low, when the amplitude is increased the dislocation 
vibrates with high amplitude in the same frequency. But the same time 
both dislocations move continually towards each other till time for 
annihilation. 

• When the dislocations move beyond the time of annihilation (in our 
case that is unit), the dislocation shows mini- vibrations except the 



6.1 Conclusion 


66 


periodicity by external stress. This vibration is due to iteration of l/2m 
of equation 3.16, which fails to converge after annihilation time. 

2. In case of 7i dislocations in a line the resultant stress field is such that one 
positive dislocation moves to +oo and one negative dislocation moves to 
— oo. remaining dislocations vibrates around their initial coordinates. 

3. In case of n dislocations on two-dimensional field they make two vertical 
clusters of of positive dislocations and negative dislocations which move with 
time away from each other.we have seen from section(2.5) that the 
dislocations are stable when ±x = 0, so they make clusters of same sign 
dislocations. 

• The vertical walls move away from each other with time. This is due to 
the fact that the net elastic energy can be reduced if they form walls 
and move away from each other. We can understand this from the fact 
the the forces acting on the dislocations from the other dislocations can 
be thought of as the derivative of the interaction energy of the two 
dislocations. 

• If the xmin is increased successively, the dislocation wall tend to 
defuse. We can understand this as, if the same no. of dislocations are 
distributed on larger field then the resultant stress field is not so intense 
to generate the regular pattern. 

• If there is no intrinsic length scale, then the walls separate out to the 
left and right. If there is indeed an intrinsic length scale then we see 
that the wall splits into the form of a sin wave. But if we make the 
length of the simulation cell in the y direction longer, the wavelength of 
the wave is also larger. What this means is that if the simulation cell 
size is much larger compared to the intrinsic length scale, then we see 
two walls to the left and right. 

• we can use this to explain the patterns observed in experiments. If the 
dislocation density is large, then the intrinsic length scale is important. 
Thus we will see the dislocations going into small pockets rather than 
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walls as we see in the simulations with a length scale. Thus, 
qualitatively what we see in the simulations is in over all agreement 
with experiments. 

4. In case of external stress the dislocation-walls are driven by external stress, 
they vibrate with time on x-axis. The dislocation- walls deteriorate with time 
also. 


6.2 Scope for Future Work 

1. In the present study, we have considered only the external cyclic stress of 
pure shear type, we can also consider the different type of loads like 
compressive , tensile and mixed state of all these. 

2. In the whole study we have considered only pure edge dislocations. We have 
scope for doing this simulation for pure screw dislocations and mixed 
dislocations also. 

3. In case of two-dimensional simulation we have neglected the inter-slip motion 
of dislocation, So we neglected the force in y direction. We can include climb 
for further study. 

4. We have not included the annihilation in our computer- code, The code can 
be modified for annihilation. 

5. If very powerful computational system is available then three-dimensional 
simulation which is more realistic, can be done. 
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